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We investigate the isocurvaton model, in which the isocurvature perturbation 
plays a role in suppressing the curvature perturbation, and large non-Gaussianity 
and gravitational waves can be produced with no isocurvature perturbation for dark 
matter. We show that in the slow roll non-interacting multi-field theory, the isocur- 
vaton mechanism can not be realized. This result can also be generalized to most of 
the studied models with generalized kinetic terms. We also study the implications 
for the curvaton model. We show that there is a combined constraint for curvaton on 
non-Gaussianity, gravitational waves and isocurvature perturbation. The technique 
used in this paper can also help to simplify some calculations in the mixed infia- 
ton and curvaton models. We also investigate possibilities to produce large negative 
non-Gaussianity and nonlocal non-Gaussianity in the curvaton model. 
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1 Introduction 



Inflation has been remarkably successful in solving some puzzles in the standard hot 
big bang cosmology [H [21 [3l H] . Inflation also predicts that fluctuations of quantum 
origin were generated and frozen to seed wrinkles in the cosmic microwave background 
(CMB) [HIE] and today's large scale structure [H [HI [HI [lOl [H] . 

In spite of the success, inflation also faces some naturalness problems. One of 
the problems is why the inflaton potential is so flat, leading to typically 10^ e-folds 
instead of 50 ~ 60, which is needed to solve the flatness and horizon problems. 

Since the invention of inflation, a great number of inflation models were pro- 
posed. Selecting the correct inflation model has become one of the key problems in 
cosmology. The currently observed quantities such as the power spectrum and the 
spectral index are not adequate to distinguish the inflation models. However, luckily 
some more quantities are expected to be measured accurately in the forthcoming ex- 
periments. For example, non-Gaussianity, isocurvature perturbation, and primordial 
gravitational waves. 

Non-Gaussianity characterizes the departure of perturbations from the Gaussian 
distribution. To characterize this departure, the non-Gaussian estimator f^i is often 
used. Using the WMAP convention, fjy^ can be written as [12] 

C = C. + IinlC, , (1) 

where ( is the curvature perturbation in the uniform density slice, and (g is the 
Gaussian part of (. This particular model of non-Gaussianity is called the local shape 
non-Gaussianity. The simplest single fleld inflation models predict that /atl < C>{1). 
So large f^L indicates a departure from these simplest models, such as the curvaton 
[T3] and (non-inflationary) ekyrotic [T?] models. To describe more general shapes, 
such as the k-type[15] and the DBI shapes |16j, one needs to calculate the 3-point 
correlation functions. 

Recently, there are hints from experiments that the local shape non-Gaussianity 
may be large. [17] claims that /atl = is excluded above 99.5% confidence level. In 
the WMAP 5-year data analysis, it is shown that the expectation value of Jnl using 
the bi-spectrum method is /nl = 51, but /nl = still lies within the 2a range. If 
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the non-Gaussianity is confirmed in the WMAP 8-year or the Planck experiment, it 
will be a very powerful tool to distinguish between inflation models. 

Another possibility beyond the simplest single field inflation model is isocurvature 
perturbation. The existence of isocurvature perturbation indicates that there are 
more than one scalar degrees of freedom during inflation. This can arise from multi- 
field inflation [181 [IH]; niodified gravity [20j, or some exotic matter content during 
inflation [21]. The WMAP 5-year -l-BAO+SN bound on isocurvature perturbation is 
a_i < 0.0037 (95% CL). So no evidence for isocurvature perturbation is shown. This 
result can be used to constrain models such as the curvaton model. 

The primordial gravitational waves also provide an important probe for the early 
universe. The amplitude of gravitational waves varies greatly in different inflation 
models. For example, chaotic inflation predicts a tensor-to-scalar ratio r ~ (9(0.1). 
While most known stringy inflation models predict r < O{10^^). The WMAP 5-year 
result, combined with BAO and SN gives r < 0.2 at 95% confidence level. This 
has put a tight constraint on chaotic inflation models. On the other hand, if future 
experiments show that r > O{10~^), it will be a challenge for string cosmology. 

One attempt to solve the flat potential problem of inflation, and to produce a large 
tensor-to-scalar ratio is proposed in [22]. The idea is to suppress the perturbation 
outside the inflationary horizon. This scenario looks like the curvaton scenario, while 
the second scalar field only creates isocurvature perturbation. So we call this scenario 
"isocurvaton" . In this paper, we show that if the isocurvaton scenario can be realized, 
large non-Gaussianity can be produced, without producing observable isocurvature 
perturbations. 

However, it is not easy to realize the isocurvaton scenario. In [22|, it is shown 
that in the simplest slow roll inflation model, with m^x^ type isocurvaton, the above 
scenario does not work. In [23] it is argued that the scenario does not work either 
for more general cases. In this paper, we prove a no-go theorem that during slow roll 
inflation, if the isocurvaton does not interact with the inflaton, then the super-horizon 
perturbation can not be suppressed. The proof also goes through for the k-type [15] 
isocurvaton, where the kinetic term of the isocurvaton and inflaton are allowed to 
be generalized. After proving the no-go theorem, we discuss some possible ways to 
bypass the theorem. 
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The proof of the no-go theorem also provides some new insights into the com- 
bined inflaton and curvaton fluctuation [24j. It is shown that in the uniform inflaton 
density slice, the curvaton propagates freely after the quantum initial condition is pro- 
vided. This provides a simplified treatment for the combined inflaton and curvaton 
fluctuation. 

We also combine the results of non-Gaussianity, isocurvature perturbation and 
gravitational waves to constrain the curvaton model. It is shown that the inequality 

\2/3 

fNL< iff^r ( ^ ) should be satisfied for the non-Gaussianity, gravitational waves 
and the temperature of the universe when CDM was created. 

This paper is organized as follows, in Section 2, we show the virtues of the isocur- 
vaton model, this is why the model is worth to investigate. In Section 3, we prove a 
no-go theorem for the isocurvaton scenario in the slow roll non-interacting multi-field 
models. In Section 4, we extend the proof to the case with generalized kinetic terms. 
In Section 5, we discuss the implications for the curvaton model. In Section 6, we 
provide a simplified analysis for the combined inflation and curvaton perturbations. 
In Section 7, we discuss some other possibilities including large negative /tvl and large 
equilateral non-Gaussianity. We conclude and discuss some possibilities to bypass the 
no-go theorem in Section 8. 



2 Virtues of Isocurvaton 

In this section, we discuss the virtues of isocurvaton. We show that some features of 
this scenario are rather attractive. That is why a no-go theorem is valuable to mark 
the forbidden regions. 

We use (f to denote the inflaton and use x to denote the isocurvaton. It is shown 
in [25] that if there is no interaction between these two components, the curvature 
perturbations for inflaton and isocurvaton on their uniform density slices are sepa- 
rately conserved. The proof is reviewed briefly in the appendix. These curvature 
perturbations can be written in a gauge invariant form as 

C, = -^-H^-^, Cx = -^-^^, (2) 
Pf Px 

where is the metric perturbation. The explicit definitions for ifj, 5p^ and 5p^ are 
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given in the appendix. The total curvature perturbation takes the form 

^ = = + r^-^. (3) 

P Px ~^ Pv 

Note that if the inflaton and the isocurvaton have different equations of state, r 
will vary with respect to time. In this case, ( is not conserved. Especially, during the 
epoch that inflaton has decayed to radiation and the isocurvaton oscillates around its 
minimum, we have oc and oc a~^. In this case, r increases with time until 
the isocurvaton decays. After curvaton decays to radiation, r is a constant, and ( is 
conserved. If we would further assume ^ (1 — r)^^ when the isocurvaton decays, 
then 

C = (1 - r)g . (4) 

From (jll) we observe that if the isocurvaton decays very late so that 1 — r ^ 1, then 
the super Hubble horizon perturbation is suppressed. 

The direct consequence of suppressing the super-horizon perturbation is to provide 
a solution to the problem of the flatness of the potential. This can make inflation 
more natural, this is because a large scalar type perturbation usually implies a non-flat 
potential. 

The isocurvaton also serves as an amplifier for non-Gaussianity, this is because if 
the initial inflaton fluctuation is larger, it should generate a larger non-Gaussianity 
than the standard scenario. This can be seen explicitly by writing 

3 ^ 

(ip = (ipg + -fNLipCpg (5) 

Combining ([T]), (jl]) and for the observable non-Gaussianity, we get 

JnL = -. fNLip ■ (6) 

1 — r 

When 1 — r ^ 1, /tvl can be large in the isocurvaton model. 

For more general shape of non-Gaussianity, the 3-point function of ip is also am- 
plified by isocurvaton. To see this, note that 

(CklCk2Ck3) = (1 - rf {Cp\,lCpV.2CvV.-i) ■ (7) 
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Generally, one can rewrite the above 3-point functions using the bi-spectrum ex- 
pression 

(Ck.Ck.Ck3) oc 6\k^ + k2 + k3)P^ViT'°'''^^(fci, h) , (8) 

where is the dimensionless power spectrum of C, and ^(/ci, A;2, A^s) describes the 
shape of the non-Gaussianity. 

So for general shape non-Gaussianity we have 

(nonlocal) 1 (■(nonlocal) /r.\ 

JnL — _ J. JNLip ■ 

As there are two components in the model, it is natural to ask whether the isocur- 
vature perturbation is produced in the model. The treatment for isocurvature per- 
turbation is the same as for the curvaton model. If dark matter is produced after 
the isocurvaton decays, the model can be consistent with the experimental results on 
isocurvature perturbation. 

In the isocurvaton scenario, the scalar perturbation is suppressed, however, the 
tensor perturbation is unaffected. This results in an enhancement for the tensor-to- 
scalar ratio. In this scenario, the observed tensor-to-scalar ratio becomes 

where tto is the tensor-to-scalar ratio without the isocurvaton dilution. If the isocur- 
vaton scenario works, and future experiments detect gravitation waves, then isocur- 
vaton with 1 — r <^ 1 can be a way to save a large number of string inflation models, 
which make the prediction that the gravitational waves are too small to be detected. 
This possibility is investigated in detail in [22] . 

However, unluckily, as we shall show in the following two sections, under reason- 
able assumptions, the isocurvaton scenario can not be realized. 



3 A No-Go Theorem for Isocurvaton 

As stated in the introduction, there are theoretical obstructions to constructing the 
isocurvaton model. In this section, we prove a no-go theorem that the isocurvaton 
model can not be realized in the standard non-interacting slow roll double field models. 
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We shall prove that in the 6p^p = gauge, to good approximation, 5x propagates 
freely, and do not feel the inflaton fluctuation or the gravitational potential. With 
this result, the gauge invariant curvature perturbation can be written as 

C, = -^, Cx = -V^-^^ = -^ + ^ + ^^5x, (11) 

Px 3x 3 

where as we shall prove, 6x is an independent stochastic source other than ip. So 
the —ip term in can not be canceled without fine-tuning , and <^ (^p can not 
be naturally realized. This result indicates that the model considered above can not 
realize the isocurvaton scenario. 

To prove the free propagation of Sx in the Sp^ = gauge, we first show that outside 
the horizon, 6x propagates freely without gravitational source term. After that, we 
show that the initial condition for 6x is determined by the quantum fluctuation before 
horizon exit, and the influence from the inflaton fluctuation and the gravitational 
potential can be neglected. 

We start with the familiar Newtonian gauge perturbation equations. Before cur- 
vaton dominates the energy density, the perturbation equations takes the form 

- 3i/(i/^/'(°) + V-^"^) - ^V'^''^ = 4vrG5pIf^ , (12) 

{2H + 3i72)^(") + 4iy^(") + = AixGdp^^^ , (13) 

+ 3i/5x^"^ + V^^5x^^^ = -2K,^(°) + 4x^('^) , (14) 

where the superscript "(n)" denotes the Newtonian gauge. The gauge transformation 
from the Newtonian gauge to the = gauge can be written as 

r (n) 

^(°)=^-i7/3, 5a;(") = 5a; + , /? = , (15) 

where x = x{t) denotes a background scalar field, and 6x stands for its perturbation. 
We assume that = Py,(p<^), so that in the 6p^ = gauge, we have 6p^ = 0. This 
assumption holds for the ideal fluid without intrinsic isocurvature perturbation, as 
well as the scalar field outside the inflationary horizon. The proof for the ideal fluid 
is straightforward, and the proof for the scalar field is given in the appendix. 

We first consider the k <^ aH limit. Changing the equations into the 6p^ = 
gauge, one can simplify Eqs. ( fT2l) and ( fT3ll as 

^ = 0, $ = iIj-H/3, (16) 
7 



note that in this paper, all perturbation variables without the superscript (n) denote 
perturbations in the 5p^p = gauge if not stated otherwise. 

Writing Eq. (JT4l) in the = gauge, and using Eq. (JT6l) . we find 

Sx + 3HSx + V^Jx = . (17) 

The i/j and (3 terms are canceled in this equation. In other words, in this gauge, x 
does not feel the gravitational potential and propagates freely. 

This result can be obtained in a simpler way. One can show that Sx is propor- 
tional to the isocurvature perturbation. Then from the well-known result in double 
field infiation that isocurvature perturbation propagates without source outside the 
horizon, we obtain that 6x propagates freely. However, we still write down the deriva- 
tion explicitly, because this derivation is rather general, holds after (f decays, and can 
be used to simplify some calculations in the curvaton model. 

Now let us consider the k ^ aH and k ~ aH case, and see whether 6x can 
feel the gravitational potential. Note that the super horizon analysis only requires 

= P<fi{Pifi), and does not require detailed information about the infiaton. While to 
investigate the horizon crossing, we need to focus on the standard single field infiaton 
plus the isocurvaton. 

We employ the results in the double field infiation model [18] to rewrite ip and x 
into the infiation direction a and the isocurvature direction s, 

X 

Sa = cosOSlp + sm95x , 5s = —sm95Lp + cos95x , sin6' = = . (18) 

Note that 6s is automatically gauge invariant. During infiation, if x ^ 0, we have 
~ during infiation. The infiation direction does not change and the isocurvature 
perturbation is obviously sourceless. However, we do not limit to this case, because 
we only require p-^ <^ during infiation. 

The perturbation equation for the isocurvature direction can be written as 

5s + 3H6s +(^ + Vss + SeA 5s = -— ^V'^"^ • (19) 

Now we shall prove that the RHS of Eq. (|T9l) is much smaller than a typical term 
in the LHS. To see this, we first estimate the fiuctuation amplitude of the infiation 
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direction. From the perturbation in the Newtonian gauge ip^"^^ + Hip^"^^ = AirGaSa^^^ 
and the slow roll condition, we have 

47rG(T25(T(°) 



< 



From the amplitude for 6a in the = gauge |18j, we have 



(5a(°) ~ (Ja^"") + —ij^''^ = {6a) 



^ r \- ' ;^=0 gauge 

Next, from the perturbation equation ip^^^ + Htp^^^ = A7rGa6a^^\ we have 



(n)| 



1.-3/2^ 
p 



(20) 



(21) 



(22) 



where we only want to count the orders in slow roll parameters, so the numerical 
coefficients are neglected. The source term in the RHS of ( |T9l) takes the form 



e e 



(n) 







- X'P 




02 


+ X' 



a 27100? if^ + ci^ 

where we have used the slow roll approximation 



(23) 



X^ - XV 



< 



XV - XV 



20X 



< 



V 



20 



02 _^ ^2 

However the quantum initial condition of 6s is \6s 
a typical term in the LHS of f|T9|) . 



+ 



X 



2X 



< H 



(24) 



~ a so when k > aH, for 



k' 



:6s 



> H%k-'l' » 



e k^ 



a 27rGa2 



(n) 



(25) 



Since the horizon exit does not take too many e-folds, the suppression in slow-roll 
parameter indicates that the initial condition of 6s is prepared by the quantum fluc- 
tuation, and the influence from the gravitational potential can be ignored. Physically, 
this result originates from the fact that the inflaton and the isocurvaton fields couple 
weakly due to the slow roll conditions. 

Note that we have ignored the back-reaction of 6s on the infiaton direction. This 
approximation can also be verified using the slow roll approximation. 

After horizon exit, in the 6p^ = gauge, ~ [25j. So in this case, 6s ~ cos 96x- 
Combining the above discussion for 6s, we conclude that the initial condition for 6x 
is prepared by its quantum fiuctuation. 
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Finally, recall (fTTI) . we get to the conclusion that 5x has its independent quantum 
initial condition, evolves freely and does not feel the gravitational potential. 

This proof of the no-go theorem can be generalized directly to the non-interacting 
multi-field isocurvaton case. So increasing the number of fields does not make things 
better. 

There are two exceptions where the above proof does not apply, namely, the 
vacuum energy and a field with completely flat potential as the isocurvaton. However, 
neither of them can serve as isocurvaton. For the vacuum energy, we can think of it 
as a shift of the infiaton potential, so it can not dilute the infiaton perturbation. For 
a field with completely fiat potential, the solutions for both a and da have a constant 
mode plus a decaying mode. The constant mode does not contribute to Co-- So up 
to the decaying mode, the flat potential case is the same as the vacuum energy case, 
and can not serve as isocurvaton. 

In this section, we have provided a direct and self-contained proof for the no-go 
theorem on isocurvaton. In order to link the double field infiaton to observations, 
analysis similar to the above proof has been performed in the literature. In a series 
of papers [26], the authors proved that the cross correlation between the adiabatic 
and entropy modes is suppressed by the slow roll parameters, and the primordial 
adiabatic mode is related to the adiabatic mode at horizon crossing by 



where T-jis is the transfer function from entropy mode to adiabatic mode. From this 
relation, one can see that the super horizon perturbation can not be suppressed in the 
context of slow roll inflation. This reasoning also extends to the interacting double 
field theory with some additional slow roll assumptions. 

4 Proof for Generalized Kinetic Terms 

In this section, we try to generalize the no-go theorem in the last section to the case 
of generalized kinetic terms. As done in the last section, we first investigate the super 
horizon evolution, and then study the horizon crossing. 

Consider the isocurvaton Lagrangian P = P{X{x), x)? where X{x) = \g^'^d^x9vXi 




sin 




1 



(26) 
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and a general dominate component originating from the inflaton p^p = Pip{pip)- In the 
k <^ aH hmit, the couphng equations for p^, and (fT2|) . (fT3|) are not changed, so 
we still have (fT6l) . 

The equation of motion for x takes the form 

Pxg^^'^.'^uX + d^Pxg^'d^x - = • (27) 

Expanding this equation to the zeroth and first order in the perturbation variables, 
we get the background and the leading order perturbation equations in the Newtonian 
gauge, 

dt{Pxx) + 3HPxX-Px = , (28) 

(29) 

Terms such as SPj^^ in Eq. ( l30il can be expanded into more explicit forms. But 
we do not need this expansion for our purpose. Note that X, Px and P^ are scalars 
under the gauge transformation (ITSl) . Using the background equation of motion (!28|) . 
it can be shown that in the (5p^ = gauge, all the source terms in (l30l) vanish, 

xSPx + PJx + xS'Px + PxSx + SHxSPx + SHPxSx -SP^ = . (30) 

Again, this result is not surprising, as isocurvature perturbation should be source- 
less after horizon crossing. 

For horizon crossing, we focus on the model that the inflaton and the isocurvaton 
have a unified generalized kinetic term. The Lagrangian of the model takes the form 

P = P{X, ^,x) , X = ^G/jV^/V'^yp-' , (31) 

where Gjj is the metric in the field space, and I, J = {1, 2} such that ip^ = ip, ip'^ = x- 
Using the results obtained in [27], the isocurvature direction perturbation (ITHj) 
can be written as 

where 

S = ^ ((1 + (^l)Ps - clPxs^^) , (33) 
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and 




1 



+ 2 



PxsPs 



Px + 2XPxx 



Px 



(34) 



where R is the scalar curvature in the field space. Note that S plays the role of 6 in 
the standard kinetic term case, which characterizes the coupling between the infiaton 
direction and the isocurvature perturbations. 

When H = 0, the isocurvaton propagates freely at horizon crossing. So we still 
have that C,^ equals C,^ plus an term originating from an independent random quantum 
initial fluctuation, and the isocurvaton scenario does not work. 

If H is large enough to provide source to the isocurvature perturbation, the above 
proof breaks down. However, the S 7^ case has not been investigated analytically 
in the literature, with only numerical results available (see, e.g. |28|). We are not 
able to prove the no-go theorem in this case. 

Note that in the proof of the horizon crossing, the generalized kinetic term will 
introduce interaction between and x- However, our proof makes no reference to the 
conserved quantities during inflation. So the interaction is not an obstruction of our 
proof. 

5 Implication for the Curvaton Models 

In this section, we flrst show that the no-go theorem proved above does not rule 
out the curvaton scenario. We also discuss some physical constraint for the curvaton 
model from the non-Gaussianity, isocurvature perturbation and gravitational waves 
experimental data. 

Let us see why the curvaton scenario is not affected by the no-go theorem. Most 
of the calculation in the above two sections applies to the curvaton scenario by setting 
'?/'*^°^ = 0. But in the curvaton scenario, it is the curvaton fleld, not the inflaton fleld, 
which produces the primordial perturbations. In the curvaton scenario, C,^p is small, 
and does not need to be canceled by the curvaton fleld. So the no-go theorem does 
no harm to the curvaton scenario. 

However, the observables we consider, namely, non-Gaussianity, isocurvature per- 
turbation and gravitational waves do put a tight constraint on the curvaton scenario. 
In the remainder of this section, we shall combine the "non-Gaussianity + gravita- 
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tional waves" [29] and the "non-Gaussianity + isocurvature perturbation" constraints 
[HI [301 [21] to provide a more complete constraint for the curvaton model. 

We first quickly review the result in [29]. Consider the simplest curvaton model. 
To distinguish it from the inflaton direction used in the above sections, we denote the 
curvaton field by x- The local shape non-Gaussianity is related to the ratio of energy 
densities when curvaton decays 

Note that by writing this equation, we have assumed f^L > 1- Otherwise, the order 
1 and order r terms in /^vl can dominate the expression. 

The curvaton starts to oscillate after inflaton decays into radiation. As the cur- 
vaton is much lighter than the inflation scale, the field value of the curvaton field is 
practically unchanged from the time of horizon exit to the time the curvaton starts 
to oscillate. So when the curvaton starts to oscillate, we have 

X 

H = m , p^ = -m'^xl , = Sm'^M^ , (36) 

where m is the curvaton mass, and x* is the curvaton field value at horizon exit. 

Another important time scale in the curvaton scenario is the time when the cur- 
vaton decays. When the curvaton decays, we have 

H = T , p^ = SM^ . (37) 

From oc and p^p oc a~*, we have = {m/V)^^'^ p^ when the curvaton 
decays. So 

.2 1/2 



6M2 V r 
In terms of /atl, we have 

2 xi \mj 
On the other hand, from the power spectrum 

pV2 = — r— = . (40) 

^ 3iT X* 12vr fNL X* 



Use (139|) and (140|) to cancel the unknown x*, we have 



/,,= 5 ey'\ (41) 



432 \rJ 
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where tt = Pt/P(; is the tensor-to-scalar ratio, and Pt = 2iJ^/(7r^Mp) is the tensor 
mode power spectrum. 

Finally, note that the curvaton decay rate F should be larger than the decay rate 
via gravitational coupling. 



m 



r>r,-^. (42) 

p 

The above results in this section have been present in [22]. Now we use them 
together with the isocurvature constraint to derive a new inequality. Use fH2]) to 
cancel m in (HTl) . we have 

fifeV'\ (43, 



432 V r , 

It is shown in [6l [30l [31] that to avoid a large isocurvature perturbation, cold dark 
matter (CDM) should be produced after the curvaton decays. So we have 

r > HcBM , (44) 

where Hcbm denotes the Hubble parameter when CDM is produced. Hcbm can be 
related with the temperature of the universe T when CDM decays as 

HcDM ^ T^Mp . (45) 

Using Eq. fj43l) . we have 

This bound should be used combined with the bound given in [29j . 

f^L < 5224/' . (47) 

The more constraining one should be used to get the final constraint. 

For example, assume /nl is produced by the curvaton scenario. If /atl = 50, then 
rr > 10"^. If the lower bound tt = 10~^ is saturated, we get T < 10""Mp ~ lO^GeV. 
This constraint is not very tight for mcDM- However, it already rules out some CDM 
candidates within the curvaton scenario, such as invisible axions, magnetic monopoles 
and pyrgons. 
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Figure 1: In this figure, we show the constraint P6|) and P7|) . The up plate of the 
figure denotes the allowed region of /atl as a function of from fHTl) . The shaded 
region is the allowed region. The down plate of the figure denotes the temperature 
when CDM decays which saturates the inequalities. In the shaded region, (1471) is 
more strict than (H6l) . In the unshaded region, (H6!) is more strict than P7|) . 

On the other hand, if Jnl = 5 then tt > 10^®. This seems more natural in the 
small field inflation models. In this case, if the lower bound is saturated, we get 
T < 10^^"^ Mp ~ lOGeV. This further rules out some dark matter candidates such as 
primordial black holes. 

The inequality (H6l) is applicable until f^i ~ 1- After that, other corrections for 
Inl begins to dominate. In this limit, > 10^^^. If this bound is saturated, then 
T < 10~^^Mp ~ 100 MeV. This bound becomes borderline for LSP and quark nuggets 
type dark matter. 

Finally, we would like to compare our result with a recent paper [31] . In [31] , the 
author also aims to get a bound for the temperature of /nl, and T. The difference 
between our work and [31] is that we use different methods to constrain the curvaton 
mass m. [31] uses the spectral index, which gives m < O.IH^. While we use Eq. (142!) 
to constrain m. This difference leads to different final results. In [31], the constraint 
takes the form 

T < 1.9 X lO-Vj/^T^iMp . (48) 
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When Trl^ < 0.023/7VL, our bound Eq. (1461) is more tight than Eq. (I48p . When 
> 0.023/7VL, Eq. (HHD is more tight. 

6 Improved Treatment for Combined Inflaton and 
Curvaton Perturbations 

Using the techniques developed in this paper, we are able to simplify some calculations 
for the mixed inflaton and curvaton scenario in the literature. 

Note that the perturbations in the curvaton scenario (setting ip^^^ = 0), isocur- 
vaton scenario, and mixed inflaton and curvaton scenario are the same. The sim- 
plifications arises because we have chosen the 6py, = gauge. This gauge choice, 
or rearranging the variables, can diagonalize the perturbation equations outside the 
horizon. 

Consider the scenario investigated in [21], with the curvaton evolving in the ra- 
diation dominated era. Using the method in [21], one is forced to solve the Bessel 
equation with source term 

S'x'''^ + SHS'x'''^ + m^6x^''^ = 4xV'^''^ - 2m\iP^'''> , H = ^. (49) 

Although this equation can be solved analytically, it saves some calculation and be 
easier to generalize if one rewrites this equation in a sourceless manner. This simpli- 
fication is just what we have done in Eq. (fT7|l . where we have considered a general 
potential for the curvaton and a general equation of state for the inflaton. 
For the m^x^ type curvaton potential, Eq. (15) reads 

5x + ^HS'x + m^Sx = . (50) 

Note that this equation has the same form with the background evolution for x- So 
without solving any differential equations, we know that for the non-decaying solution, 

Sx=—^Xo, (51) 
Xo 

where 6xo is an integration constant. Use (HM to translate this result into the New- 
tonian gauge, we have 

Sx^""^ = - (Sx'o'^ - Xotoi^^^^) + . (52) 
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Note that 6x0 — Xo^o'^o J^^^ ^ constant, which can be redefined to Sxo ■ Then 
we recover one of the key results in 

= -Sx^o'^ + , (53) 

Xo 

where we have not made any reference to the equation of state of the infiaton com- 
ponent. Other results for the m^x^ potential in [24j can be recovered similarly. 

7 Generalizations for the Curvaton Model 

In this section, we investigate the possibility for large negative fj^ji and nonlocal shape 
fisiL in the curvaton model. This possibility can be realized by phantom curvaton 
and k-curvaton respectively. These models seem exotic, and not supported by any 
evidence so far. For example, the vacuum stability and the quantization problem in 
the phantom model are not solved. However, we investigate these models as pure 
phenomenological possibilities. 

In reference [6], the WMAP5 data has been analyzed by two different methods. 
It is intriguing that the two methods prefer central values of /^£"' opposite in sign. 
In particular, both at the 95% confidence level, the bispectrum analysis gives the 
best estimate —9 < /j^£"' < 111, while the analysis of Minkowski functionals prefers 
— 178 < /^£"' < 64 in contrast. It is still unclear why they are so different. However, 
if one naively disregards the bispectrum analysis for the moment, and takes seriously 
the central value from Minkowski functionals, we will be motivated to search for 
models with ^ — 1. Let us see what will happen if the curvaton component is 
phantom-like [32]. In that case, the Eqs. ([3]) and ([2]) are still valid. However, now 
p^> 0. So we have r < 0, and has the different sign from (. It is well-known that 
a different sign in ( should produce a different sign in /^vl- The calculation for f^^ 
goes through in the phantom model, so when |r| <^ 1, we have 

/iVL ~ - < -1 . (54) 
r 

Although the bispectrum analysis of WMAP5 (which is widely taken as the best 
estimate) does not prefer a large negative fjyi, its lower bound /^^ ~ — 9 still allows 
for the phantomlike curvaton to live in a narrow space. From the opposite viewpoint, 
this can be another piece of evidence that nature disfavors phantom. 
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It is also worth to note that if the index of equation of state for the curvaton crosses 
— 1 [33], then the non-Gaussianity produced by the curvaton model also crosses — 1. 
To realize this possibility, one usually need more than one curvaton fields |34j . 

Now consider the k-curvaton possibility. If the curvaton has generalized kinetic 
terms, then the equilateral non-Gaussianity for the curvaton is also large. Similar to 
([9]), we have 

(nonlocal) 1 (nonlocal) {e.t^\ 

Jnl — ~Jnlx ■ y^"^) 
This amplification can easily produce very large equilateral non-Gaussianity. Note 
that /^^'^^°'^^'^ ~ For example, if l/c^ ~ 5, and /tvl — 50, then we find 

^^noniocai) _ 250. The experimental bound -151 < /^T^"""'^ < 253 (95% CL). If both 
large local and nonlocal /^^°'°'^'^^'' is observed, the k-curvaton provides a satisfying 
explanation. 



8 Conclusion and Discussion 

To conclude, in this paper, we have investigated the isocurvaton scenario. We found 
that although the isocurvaton scenario possesses attractive features such as enhance- 
ment of non-Gaussianity and gravitational waves, the scenario can not be realized in 
the slow roll multi-field models. This no-go theorem can be extended to generalized 
kinetic terms with assumption S = 0. The techniques used in this paper can simplify 
some calculations in the mixed curvaton and inflaton scenario, providing an easier 
investigation for more general mixed perturbations. 

We showed that the no-go result does no harm to the curvaton scenario. However, 
the experimental bound on non-Gaussianity, isocurvature perturbation, and gravita- 
tional waves provide a combined constraint (Eq. (H6l) ) on the curvaton model. 

We also investigated the phenomenology of phantom and kinetic curvatons. We 
showed that the phantom curvaton provides fNL ^ —1, and the k-curvaton provides 
very large equilateral non-Gaussianity as well as the local non-Gaussianity. 

Finally, let us discuss some possibilities to bypass the no-go theorem for isocurva- 
ton. The following possibilities are not covered by the no-go theorem: 

1. Adding interactions. It is reported from numerical calculation that interaction 
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can suppress the super horizon perturbations [35]. It would be interesting to 
investigate whether similar mechanisms can realize the isocurvaton scenario. 

2. Relaxing the slow roll condition for the isocurvaton field. It is challenging to 
construct fast rolling isocurvaton field which can fit the experimental results. 

3. Other form of generalized kinetic terms, including separately generalized kinetic 
terms for infiaton and curvaton during infiation, high derivatives like the box 
term. The possibility S 7^ is also worth investigating. 
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Appendix 

We explain the notation we use, and review some well-known facts in the cosmological 
perturbation theory. 

In the linear perturbation theory, assuming a fiat universe {K = 0), and without 
choosing any gauge, the metric for the scalar perturbation takes the form 

^'^'^ \ -/3,^ -a\{l-2^)5^, + 2E,,,) 
The Newtonian gauge is defined by setting 

For the Sp^ = gauge, the equation Sp^ = is just one gauge condition, and as 
in [36], we set the other gauge condition to be = 0. In this notation, the gauge 
transformation takes the form of Eq. ( fTSll . 

The conserved quantity can be introduced as follows. Assuming that there is no 
energy change between ip and the local energy conservation equation for ip takes 
the form 

i7-^ = -i 9,{p, + 6p,) 

3 p^ + Sp^ +Pip + Sp^ 
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Note that if is a background quantity, and does not change with the spatial coordi- 
nate. If we assume the pressure is a function of only the energy density, then in the 
5p^p gauge, we have = 0. Thus the RHS of (l58ll is also independent of spacial co- 
ordinates. In order that (1581) holds, ip should also independent of spacial coordinates 
outside the horizon. 

As a perturbation variable, ip should have no zero mode, so does ip. The only 
possibility is = 0. We can define a conserved quantity 

C. = -^Lp,=o' (59) 

which is conserved after horizon crossing. This conserved quantity can be rewritten 
in the gauge invariant form 

= - if% . (60) 



From the same reasoning, there is also a gauge invariant conserved quantity for 

C^ = -^-H^-^. (61) 
X 

Note that these conserved quantities can be defined beyond the leading order 
perturbation theory. But we only need the leading order result for our purpose. 

The above proof is under the assumption that in the 6p^ = gauge, we have 
Sp^ = 0. This assumption is obviously true for fluids such as radiation and matter. 
It is also worth to note that this assumption is also true for inflaton after horizon 
crossing. This is because the above statement can be rewritten as the adiabatic 
condition 

P^Sp^ = p^6p^ . (62) 

This condition can be checked directly using the Einstein equations in the k <^ aH 
limit. 
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